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1. Introduction 

Effects of randomness have supplied fundamental 
problems in condensed matter physics and localization 
due to interference of quantum mechanical electrons are 
well studied as the Anderson localization. Although we 
have well established understanding of the localization of 
non-interacting electrons, information of the correlated 
electrons with randomness is still missing. It was mainly 
due to lack of reliable numerical techniques for the cor- 
related electrons. 

Today, for the one dimensional correlated systems 
without randomness, lots of numerical results are col- 
lected by the Density Matrix Renormalization Group 
(DMRG) 1 ) method and consistent understanding with 
analytical predictions has been achieved. In this paper, 
we plan to apply DMRG for the random electron systems 
by calculating direct responses of the system with elec- 
tric held. 2 ) At first, random systems without interaction 
are carefully investigated. Then we try to treat both of 
interaction and randomness in one dimensional systems. 



■ 2. Model and Dielectric response 



We investigate a one-dimensional spinless fermion 
model with nearest neighbor electron-electron interac- 
tion in the presence of random potential. The Hamilto- 
nian is given as 
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H = -t ^(cjc l+ i + c} +1 Cj) + V ^ nitii+i + ^ eic\ci, 

i—l i—1 i—1 

(i) 

where rii = c\ci and is random potential, which dis- 
tributes over the interval [— W/2, W/2] uniformly. We 
set t = 1 for simplicity and consider the half-filling case 
and impose an open boundary condition. In the absence 
of disorder, the system is metallic for — 2 < V < 2. But 
for half-filling case, at V = 2 the system undergoes a 
metal-insulator transition and for V > 2 the system has 
a finite charge gap. In attractive interaction region at 
V — —2, the system becomes unstable due to phase sep- 
aration. When the randomness is present, the system is 
always insulating due to the Anderson localization with- 
out interaction. The interplay between the randomness 
and interaction can be interesting and possible metallic 
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phase is expected for some range of negative V^. 3 " 5 ^ 

In this parer, we focus on the dielectric response of the 
system. In order to observe the dielectric response, we 
apply the electric held E to the system. As the second 
quantized form of the potential, —Ex, the coupling term 
He is added to the Hamiltonian. Then full Hamiltonian 
of the system is given by Ht = H + He where 



Hp 



-E 



(2) 
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As a function of E, the polarization P of the system 
is defined by 
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where Eq is a ground state energy and {ti^e represents 
the ground state expectation value of n^. 2 ^ Here we used 
the Feynman's theorem to derive the last equation. 

For a finite value of the electric field, which is compa- 
rable to the Mott gap, we expect a collapse of the local 
charge gap due to the interaction. Also in the Ander- 
son insulator, reconstruction of the charge by transfer- 
ring electrons above the tunneling barrier could occur. 
In each case, we can obtain information of the charge 
degree of freedom above its ground state. 

We turn next to the linear response regime. In this 
regime, we calculate zero-field dielectric susceptibility as 
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From the susceptibility x, we directly obtain informa- 
tion whether the ground state is metallic or an insulator. 
In thermodynamic limit, \ is diverging if the system is 
metallic, but converging to a finite value if it is an insu- 
lator. Indeed x ~ L 2 is expected by the perturbation for 
the pure non-interacting system(VF = 0, V = 0). 

In order to calculate the charge distribution and the 
ground state energy, we use the exact diagonalization 
for V = 0, and DMRG for finite V. For the application 
of DMRG, we use the extended infinite-size algorithm 
by Hida, 6 ) which enables us to treat also non-uniform 
lattice models. We perform three or four finite lattice 
sweep for the convergence. The retained states for the 
block is 60-100 to keep the truncation error to be less 
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than 10" 8 . 

3. Results and Discussion 

At first using exact diagonalization, we calculate the 
polarization P as a function of the applied electric field E 
in the absence of electron-electron interaction. The typi- 
cal P — E curve is shown in Fig. 1 and Fig. 2. The polar- 
ization P is a smooth function of E without randomness. 
As we increase the electric field from zero, P is almost 
linearly increasing and approaches to a finite value, at 
which all electrons are collected to one side. Then at the 
critical field E s ~ 1/L, P saturates as P ~ L. 

In Fig. 2, we plot the P—E curve for finite randomness 
strength W . In the presence of randomness, the P — E 
curve exhibits a stepwise behavior. This step is caused 
by crossing of the one particle energy levels. Namely, 
at some critical E, level crossing between the highest 
occupied state and the lowest unoccupied state occurs. 
By this process, the charge reconstruction of the ground 
state occurs which corresponds to the electron tunneling 
between localized states. 



for W — 3, x is convergent to a finite value which is 
consistent with the insulating ground state. For inter- 
mediate W, for example W = 0.5, saturation of x is n °t 
observed up to L = 500. The localization length £ seems 
to be larger than the system length we used. Then fi- 
nite size scaling analysis is necessary. The detail can be 
published elsewhere. 
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Fig. 3. Average logx as a function of log 1/L. The system sizes 
are between 10 and 500. (V = 0. W = 0, 0.5, 1.5, 3.0 and 5.0.) 




Fig. 1. 
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Polarization P as a function of applied electric field E. 
■■ and L = 20, 40, 60 and 80.) 




Using DMRG, we study also the systems with elec- 
tron correlation. In Fig. 4, logx versus log 1/L for 
V = —1.4 is plotted. We averaged over 64 realizations 
of the disorder potential. When the randomness is suffi- 
ciently strong, we observe the saturation, which implies 
the localized ground state. However, we need more ex- 
tensive analysis to obtain conclusive results. The detail 
analysis with finite size scaling will be given elsewhere. 
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Fig. 4. Average logx versus log 1/L for V = —1.4. The system 
sizes are 10, 12, 16, 20, 26, 30, 34, 38, 40, 44, 50 and 60. (W = 
0, 1.5, 3.0, 4.0 and 5.0.) 
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Fig. 2. Polarization P as a function of applied electric field E. 
(W = 5.0 and L = 20.) 



Next let us focus on the linear response region. In 
Fig. 3, we show the susceptibility x as a function of 1/L 
by numerically differentiating the ground state energy. 
We averaged over 500 realizations of the disorder poten- 
tial and we applied the electric field E = ±10~ 5 . For 
W = 0, x increases as x ~ L 2 and x seems to diverge. 
This means the system is metallic. On the other hand, 
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